DAMTP 94-17 



Classical Functional Bethe Ansatz for SL{N): 
Separation of Variables for the Magnetic Chain 

D.R.D. ScottQ 

Department of Applied Mathematics and Theoretical Physics 
Tlj" . University of Cambridge 

; Silver Street, Cambridge CBS 9EW, U.K. 

Abstract 

' The Functional Bethe Ansatz (FBA) proposed by Sklyanin is a method which gives separation 



o 



variables for systems for which an i?-matrix is known. Previously the FBA was only known for SL{2) 
\ and SL{2>) (and associated) i?-matrices. In this paper I advance Sklyanin's program by giving the 
■ FBA for certain systems with SL{N) i?-matrices. This is achieved by constructing rational functions 
A{u) and B{u) of the matrix elements of T{u), so that, in the generic case, the zeros Xi of B{u) are 
', the separation coordinates and the Pi = A{xi) provide their conjugate momenta. The method is 
illustrated with the magnetic chain and the Gaudin model, and its wider applicability is discussed. 
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1 Introduction 



1.1 Separation of variables in Integrable Systems 

A classical system with a 2n-dimensional phase space is said to be integrahle if there exist n inde- 
pendent functions on the phase space (called integrals of motion) which Poisson commute among 
themselves and with the Hamiltonian. 

Classically a system is said to be separable if there exist variables in which the Hamilton-Jacobi 
equation separates. Traditionally this means introducing n separation constants so that the Hamilton- 
Jacobi equation (involving all n pairs of conjugate variables) can be replaced by n equations involving 
one each (and possibly a further equation constraining the allowed separation constants at each en- 
ergy). 

Liouville's Theorem says that for integrable systems there exist action-angle variables (at least 
locally) in which the Hamiltonian is a function of the action variables alone and hence separates. 
Thus for integrable systems it is known that separation variables exist (at least locally). However the 
formula for constructing action-angle variables (which involves integrating over invariant tori) is not 
tractable in general. 

What is needed is an explicit method for constructing separation variables for arbitrary integrable 
systems. Two such methods exist. One uses algebraic geometry and has been applied to loop algebras 
with linear Poisson brackets |||]. The other, the Functional Bethe Ansatz (FBA) works for systems 
with SL{2) or SL(3) i?-matrices of certain forms. In this paper I extend the FBA to SL(N) for 

> 3. I work with the magnetic chain and Gaudin magnet, but the method and results are much 
more widely applicable as I will discuss at the end of the paper. Although this extension has equations 
in common with the algebraic geometric technique no knowledge of algebraic geometry is required to 
understand this paper. 

Some insight may be gained from considering the reverse problem (i.e. how to construct integrals of 
motion given a separation) which is much better understood. Jacobi's Theorem can be interpreted 
as saying that if a system separates then it is integrable, with the integrals of motion corresponding 
to the separation constants. Thus given a separation of a Hamiltonian there is a systematic method 
for constructing the integrals of motion, and the separation of the original Hamiltonian induces a 
separation of these integrals of motion. Thus it is natural in integrable systems to think of separating 
the system without any reference to a particular integral of motion. Separation of variables may 
then be (re)defined as seeking variables Xi and pi, with canonical Poisson brackets, which satisfy n 
separated equations, 



where in these equations the dependence on separation constants has been replaced with dependence 
on the integrals (which are of course constants of the motion). This is the starting point used by 




(1) 



Sklyanin for his Functional Bethe Ansatz. 
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1.2 /^-matrices and Lax Pairs 



The two methods aheady proposed and this paper ah use the matrix formahsm associated with R- 
matrices and Lax Pairs. To understand this paper only the basics of this formahsm need be known, 
further details may be found for example in [Q, ^, ^] . 

In the i2-matrix approach the operator content of the theory is contained in an x matrix 
function on the phase space T(u) depending on a spectral parameter u, the algebraic structure of the 
theory is given by the so called i?-matrix algebra. This can be quadratic, 

{Tiu),T{v)} = [R{u-v),T{u)T{v)] (2) 

or linear, 

{T {u),T (v)} = [R{u-v),T {u)+T (v)] . (3) 

Here T denotes T{u) ® 1 etc. and R is an A^^ x A^^ matrix acting in the tensor product. 

The i2-matrix must obey a consistency condition, the Classical Yang-Baxter equation. In the case 
of Lax Pairs we have a pair of matrices L(u) and M{u) depending on spectral parameter n, L{u) is a 
function on the phase space, while M{u) is usually a constant. L and M give a Lax representation of 
a system if the equations of motion of the system are equivalent to the following evolution equation 
for L, 

^ = [L(n),M(n)] . (4) 

From this equation it is clear that the quantities generated by the spectral invariants of L[u) are 
constants of the motion, but a priori they do not Poisson commute. The necessary and sufficient 
condition that the eigenvalues Poisson commute is that the Lax matrix obeys an i2-matrix type 
Poisson brackets Q, which for antisymmetric i?-matrices reduces to equation (^). 

This paper uses the i?-matrix formalism. What is important to know for this paper is that in 
general the spectral invariants of T{u) provide the integrals of motion. So long as they give rise to 
enough the system is integrable. 

1.3 The Method 

The search for separation variables can be thought of as happening in two stages, (1) look for vari- 
ables which give rise to separated equations, then (2) check that these variables have the correct 
commutation relations. 

Crucial to both methods is the idea that if Ci is an eigenvalue of T{ui) then Ui and Q automatically 
obey the following 'separated' equation, 

det (0 - T(n,)) = . (5) 

This is indeed a separated equation in that it only depends on itj, C,i and the spectral invariants of 
T{u) (i.e. the integrals of motion). Any such ui and Q found in this way therefore complete step (1). 
The problem is to choose them in such a way that they have the correct Poisson brackets and form a 
basis of the phase space. 
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1.3.1 Functional Bethe Ansatz 

The Functional Bethe Ansatz was proposed by Sklyanin |8|, 10 1 as a method of finding separation 
variables for integrable systems, this paper is only concerned with the classical version. It is a blending 
of ideas from the Bethe Ansatz and separation of variables. It is applicable to systems for which an 
i?-matrix (or indeed Lax) representation of a certain form is known. 

The FBA was first proposed for the case of 2 x 2 matrices. In this case the basic objects are 



C(u) D{u) 

and these obey i?-matrix relations. The Functional Bethe Ansatz says that if the {xi,i = 1, . . . , M} 
are the zeros of B{u) then the Xi provided a set of separation coordinates and the set of Pi given by 
Pi = A{xi) are their conjugate momenta (or possibly a function of them). The separated equations 
are then, 

det (Pi - T(xi)) = . (7) 

It has been successfully applied in this (2 x 2) case to a wide class of systems. The FBA for 5'L(3) 
was studied by Sklyanin in for the 5L(3)-magnetic chain and the Gaudin model. In this paper he 
conjectures how the method might be extended to SL{N). For the magnetic chain the conjecture is 
as follows (with a similar one for the Gaudin magnet but with canonical Poisson brackets). 

Conjecture 1 There exist functions A and B on GL{N) such that the following two assertions are 
true. Firstly that A{T) is an algebraic function and B(T) is a polynomial of degree D = MN{N — \)/2 
in the matrix elements Tmn- Secondly that the variables Xj,Pj (j = 1,...,D) defined from the 
equations, 

B{T{xi)) = 0, Pj = A{T{xj)) (8) 

have Poisson brackets 

{xj,xk} = {Pj,Pk} = 0, {Pj,xk} = Pj6jk (9) 
and, besides, are bound to the Hamiltonians (integrals of motion) by the relations 

det(P, -r(xi)) = . (10) 

The scheme proposed by Sklyanin to obtain A and B is to put T into block triangular form by using 
a similarity transform 

T' = K-^TK (11) 

where is a matrix depending on some parameters /ci, . . . kq. (I shall show that Q need only be A^— 2.) 
The ki should then be eliminated from the resulting equations to leave a single equation B{u) = 0. 
No similarity transformation is required for N = 2 and the similarity transformation for = 3 was 
performed by Sklyanin in Ref[^, where he proves the commutation relations for the magnetic chain 
and Gaudin magnet. 
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1.3.2 Algebraic Geometric Method 

M.R. Adams, J.Harnard and J.Hurtubise Q have used algebraic geometric metliods to find separation 
coordinates for systems that have Lax representations with linear i2-matrices and Lax matrices of the 
form, 

M J 

L{u) = uY + uY '— (12) 

with Y G gl{n) and Oj as complex constants (it also works for multiple poles). 

In this approach separation variables are constructed as the generically distinct finite solutions of 
the equation, 

M(n,C)Fo = (13) 

where M is the classical adjoint of 

M{uX)-=T{u)-CI (14) 

The solutions of this equation give Darboux coordinates in systems with linear Poisson brackets of a 
certain from. The solutions of this equation are again bound to the integrals of motion by equation (|5|). 
As a defining equation for the separation variables equation (|13|) has the advantage that (because it 
can be treated in the language of algebraic geometry) degenerate cases can be handled easily. Perhaps, 
most importantly, the commutation relations can be straightforwardly calculated. 

A disadvantage as compared with Sklyanin's scheme is that the defining equation involves two 
parameters (the spectral parameter and an eigenvalue parameter). 

1.4 Overview of Paper 

In this paper I advance Sklyanin's program by giving the FBA for SL{N). This is done by constructing 
polynomials A and B in the matrix elements, so that generically the zeros Xj of B{u) give the separation 
coordinates and the Pi = A{xi) provide their conjugate momenta in the cases of the magnetic chain 
and Gaudin model. This paper deals primarily with the magnetic chain which is reviewed in section (|2|) 
(the Gaudin model being dealt with at the end (section (Hj)).) 

In section (^) I obtain candidates for A and B for the case of A x matrices, in the sense that 
produce variables that give rise to separated equations, without however any consideration of whether 
they have the correct commutation relations at this stage. Sklyanin's program of similarity transforms 
is used. 

In section (0) it is shown that generically A and B give the same separation coordinates as obtained 



from equation (13), and then this equivalence is used to prove the commutation relations for the 
magnetic chain. In doing this the crucial calculation in Ref. 2 is extended to show that equation ( |T3|) 
gives separation variables in the case of quadratic Poisson brackets (given by the permutation R- 
matrix, see equation (pO|)). 

In section (^) it is shown how the results carry over to the (linear i?-matrix) case of the Gaudin 
Model, and to systems with more general i?-matrices. 
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2 The Magnetic Chain 

The variables of the non-homogeneous classical SL{N) magnetic chain are a, /? = 1, . . . , N;m = 
1, . . . ,M (where M is the length of the spin chain). The variables are not completely independent but 
a related by J2a=i '^^a^ — 0- They obey the following Poisson brackets 

i^Lili^ ^i2l32^ ^ ('^ii/32'^a2/3i - 'S'^2/3i'^"lft) ^mn (15) 

which define the Kirillov-Kostant Poissonian structure on the direct product of M coadjoint orbits of 
SL{N). The center of the algebra is generated by the eigenvalues A™ of the matrices S^^^ 

N N 

det - a) = n (AS^ - A) . EAa=0 (16) 

a=l a=l 

I shall fix the orbit, by taking A™ to be fixed numbers. Furthermore I shall assume that I have a 
generic orbit by requiring the eigenvalues of 5^*"^ are distinct. The manifold defined by equations ( |T6|) 
and having dimension MN{N — 1), is then equipped with a non-degenerate Poisson bracket (p!5[). 
The monodromy matrix may be defined as, 

{u- 82 + 5(2))(n-5i + 5«) (17) 

where Z is an x number matrix with distinct eigenvalues and 5m are some fixed numbers, and 
u is the spectral parameter. 

The matrix elements of T are polynomial in u of degree M (length of the magnetic chain). The 
T(n) obey the following quadratic R-matrix relations, 

{ra,/3i(n),r„2/32(i;)} = --^(r„2/3^(u)r«^/32(t;) -Tc,^fi^{u)Tc,.,(3^{v)). (18) 
or in the formalism used earlier, 

{r(n),T(f)} = [R{u-v),T {u)T [v]] 

Riu) = - (19) 
u 

where V is the permutation matrix in the tensor product i.e., 

1 ^ 

R{u) = - E tij tji (20) 

where the Cij form the usual basis for N x N matrices with a 1 in the ijth position and zeros else- 
where. The spectral invariants tu{u) of the matrix T{u) may be defined as the elementary symmetric 
polynomials of its eigenvalues, tu{u) = tr A'^ T{u),i' = 1,...,N. t]\f{u) = det(T(n)) contains the 
central elements and is therefore taken to be a constant function of u. The non-leading coefficients 
of the remaining invariants provide a commuting family of MN{N — l)/2 independent Hamiltonians 
(see e.g. pi), since this is half the dimension of the phase space the system is integrable. 
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3 A,B in the SL{N) case 



In this section I use the similarity transformation equations to show that the following are candidates 
for A and B, i.e. that the Xi and Pi derived from these equations allow separation (the proof that 
these variables obey the correct commutation relations is postponed to the next section.) 



Ainu)) 

B{T{u)) 



-ni2---jjv-i-'iiA^-'i2Af ' ' ' ^~ 



N 



i2N 



det(M) 

rpN-2 rp(N-l) 



(21) 
(22) 



In the case the case of the magnetic chain, where T{u) is a polynomial in u of degree M, B{u) is 
generically degree MN{N — l)/2 insuring that it defines the correct number of separation variables. 

In this section indices labeled i,j,k range from 1 to — 1, indices labelled m,n range from 2 to 
N — 2 and repeated indices are summed over their appropriate ranges. 

Let 



/I k2 
1 



K 







The requirement that, 



Vo 



KTK-^ 



kN-i 0\ 
... 





1/ 



(23) 



fAiu) 

# 

V # 



# 



(24) 



where # denotes arbitrary matrix elements, implies that, 

A{T{u)) = hTii 

and that the fcj's must satisfy the following set of — 1 equations, 

kiTim - kmkTu = fm = 2 - 1) 



kmkiTil 

kiTijsr 



(m 




(25) 



(26) 
(27) 



where i is summed from 1 to — 1. These are A^ — 1 equations in A^ — 2 unknowns our aim is to 
eliminate the unknowns and obtain a single consistency equation B{T{u))=0. The first N — 2 equations 
are quadratic and only the last is linear, however an equivalent set of linear equations can be obtained 
from these. Multiplying the equations ( [2^ ) by T^n and summing from 2 to A^ — 1, one obtains, 

(28) 



kjTjYnTraN k^T^j'q kjTii — 

and the last linear equation may be used to replace the last term by kiTnTi^, thus 

kjTjkTkN = 0. 



(29) 
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Now we have two linear equations, further equations can now be obtained by iterating this procedure 
(e.g. to obtain the next equation in the series multiply each of the original quadratic equations by 
TmiTiN and sum them using ([29| ) to eliminate the quadratic terms). Introducing the notation 

^=T^^;^'\ ^^^=T^k (30) 

the linear equations read, 

A:,SS = 0. (31) 

The first A^ — 2 of these equations are sufficient to determine the fej's, and may be written in matrix 
form as, 

where M is the {N - 2) x {N - 2) matrix, 

Mmn = 4^'^- (33) 

Hence, 

det(M)fc„ = (34) 

where N is the transpose of the matrix of cofactors. Multiplying the L.H.S. of equation(^) (with 
a = — 1) by det(M) and substituting, 

det(M)A:,S;^ = det(M)sS^-') + iV^^sS';;)^^-') (35) 

this is chosen to be B. Likewise A is obtained by eliminating the A^j's, giving, 

A{T{u)) = eiii2...i^_,i;^^^^]yS|^]y • • • ^l^-aA^ (^'') 
B{T{u)) = ei^i^...ij^_Jl\^j^T\^j^ ■ ■ ■t!~-^_J^t!i^_^Ij (38) 

In these expressions each index is summed from 1 to — 1 however A and B remain unchanged 
if the sums are extends to sums from 1 to A^, thus 

A{T{u)) = e.,.,...,_,T,,^i;2^...7;^-2^-^^ (39) 



BiT{u)) = 6,,,,..,^„,r,,rf^^...7;^:2 r/^_^;^ (4o) 



det(M) 

rpiN-l) 

To see this notice that 

= {T')iN + (T^-^)iivF(') + • • • + TiATF^^-i) (41) 
and that in A and B these appear within antisymmetric sums so only the first terms contribute. 
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4 Proof of commutation relations 



Summation convention is NOT used in this section. 

In this section I prove that the Xi and Pi given by B{xi) = and Pi = A{xi) have the fohowing 
commutation relations, 

{xi,x,} = {Pi,P,} = (42) 
{xi,Pj} = Pi6ij . (43) 

The proof works in two stages first I show that the {xi,Pi) can be equivalently defined as the solutions 
of another equation. Then this equation is used to calculate the Poisson brackets. 

4.1 Equivalent defining equation 

In this subsection I show that (if the B has the correct number of zeros then) then (xi, Pi) defined by 
A and B can be equivalently defined as the generically distinct finite solutions of the equation, 

M{u,C)Vo = (44) 

where M is the classical adjoint of 

M{u,C):=T{u)-CI (45) 

It is very natural to think of the separation variables as being defined by these equations since they 
do give eigenvalues and moreover define the variables in terms of things for which we can calculate 
Poisson brackets. The solutions of this equation were shown to give Darboux coordinates in systems 
with linear Poisson brackets [|[. 
To see this consider, 



M^{C) = Kiui)Miui,C)K-\ui) 



(Pi-C ••• o\ 

# # 



(46) 



V # #/ 

where # denotes arbitrary matrix elements. Factors of {Pi — C) may be pulled out of all but one 
of the non-zero cofactors, 

M\{0={* /f^"9'' ) (47) 



V {Pi- C)B,j , 

where q and Bij are the relevant (determinant) factors that multiply (Pj — (") in each cofactor. 
Thus for generic Vq (i.e. one for which Vq = 0), Zi = Pi solves the equation, 

M'{u^,Zi)Vo = (48) 

and hence X = Ui and C, = Pi are generic solutions of the equivalent equation 

M(A,C)Fo = (49) 

In the case that B does not give as many solutions as (^) this proof breaks down. Indeed if B 
does not give the same number of zeros as half the dimension of the phase space then the method can 
not be applied in its present form. For example the i?-matrices of Kuznetsov |Tl| (for the reducible 
systems of Kalnins et al. [^, |l3|) cannot be used to define separation variables in this manner. It 
would be very satisfying to have a FBA type method that could handle such degenerate i?-matrices. 
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4.2 Proof of commutation relations 

bmce are given as the generic solutions of equation ( p|) we may choose (Vo)i = diN^ where 

A^'o is fixed but arbitrary (and could be chosen to be 2, 3, . . . , — 1 or N). Then (xj, Pi) are given by 
the conditions, 

MfcjVo(A,C)=0 (50) 
Generically these points are specified by just two of these, choose 

MiNo = M2N0 = (51) 



and the matrix 



du dC 

is invertible. Hence 

{Cu,Uf,} {Cu,Cii} 

p -1 f WlNoiUiyXu), MiN^iUf,, Cfi)} {MlNo{u,y,Cly),M2No{u,^,C^J,)}\ ( pT^^^l /coN 
" I {M2No{Uu, Cu), MlNoiUf,, C^)} {M2N,{U„ M2N,{u^, Cm)} / ^ ^ ^ ^ ^ 

The Poisson brackets of the adjoint can be calculated from the Poisson brackets of M by using the 
derivation property, 

{*,(„, 0, M.(.. ,)} = E (M„(„, 0, M„.(„, (54) 

Now ^ 

{Mpq{uX),Mrs{v,r])} = {Trq{u)Tps{v) - Tps{u)Trq{v)) (55) 

U — V 

and it follows from M{u, C)M{u, C) = det (M(n, C)) / that 

dMij{u,C) _ Mgp{u,C)Mij{u,C)-M.,p{u,C)Mgj{u,C) 



dMpq{u,C) det M(u,C) 

Substituting these into ( |54|) one obtains, 

{MiNo , MkNo} = E (-^ 

qsMiNoiu) - 6isMqNo{u)j (SgqMkNoiv) - SkqMgNoi'^ 

pqrs 

- (6rpMiNo{u) - 6rNoMip{u)j [SprMkNoiv) - SpNgMkriv) 

dMiN,{u) 



(56) 



+ (SsqMkNoiv) - 6qkMsNo{v)) " Srq (^6rsMkNo{v) " 6pNoMkr{v)) j (57) 

+ (^6rq (SqsMiNoiu) - 6isMqNo{u)j - 6ps (^drpMiNoiu) - 5rNoMip{u) 

Thus by taking the appropriate limit with u ^ fi the left hand side vanishes. (After noticing that 
MaNoi'^i: Ci) = the two remaining terms vanish since g^"^" = 0.) 

{MlNo{Uu,Cu),MiNo{v^„r]^M)} = {MiNo{u,„Ciy),M2No{Vp,'np)} 

= {M2No{u^,Cu),M2No{v^^,V^^)} = (58) 
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Hence 

{u^,u^} = {C Cm} = {un, Cm) = ^ ^ ^ (59) 
It remains to calculate the Poisson bracket {uiy,(^i^}. One may use, 

{Mpg{uX),Mrs{uX)} = ^Tps - ^Tr, (60) 



Thus 



/A> f n n7 ( r\\ \^ ( 9MiNodM2No dM2No dMiNo \ dMps ^^^ 



pqrs 

and from equation (^6|) we see, 



\^ dMpg dMrs dMpg dMrs J (det M)2 \^ V ^'^ 

+ MsrMgNo {MiNoM2p " M2iVoMip) 

+ MgNoMsNo {MlpM2r - M2pMir) ) (62) 



The last term in brackets contains three zeros and never gives a contribution. The other terms also 
vanish when they are multiplied by ^^^M^q and summed over (M^g may always be multiplied by an 
M that is not of the form thereby cancelling a det M in the denominator but leaving two zeros 

in the numerator). Thus the only contribution comes from the product of the first two terms of (62) 



and C~^f^^rg- However this is equal to (^detFu. To see this recall that, 



^ . 17 f ^^1^0 ^^2jVo dM2No dMiNo \ dMpr ^ 

deti^, - 2^ I dMpg dMrs du ^ ^ 



pgr 

Thus explicit calculation shows that, 

f dMiNo dM2No dM2No dMiNo \ ( dMps _ dMp^ 



\ dMpg dMrs dMpg dMrs j \ du du 
1 dM, 



5 (detM)2 du [^^NoMvNo {MsrM2p - MspM2r) + ^N^MrNo {MspMlr - MsrMlp 

+MiN,MsNo (MrrM2p - M.^Msr) + MsiVo^.^Q (MrpMlr " MlpMrr) j = (64) 

where the terms in brackets give the required third zeros because they are the subdeterminants of 
a matrix of rank 1. Hence 

{MiNo {u, C), M2N0 {u, C)} = C det F, (65) 



thus substituting in equation ( |53| ) we find 

{uuXu} = Cu (66) 

as required. 
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5 The Gaudin Model and other systems 
5.1 The Gaudin Model 

In the Gaudin model (which may be considered a degenerate case of the magnetic chain) T{u) has the 
fohowing form, 

m=l 

where S"^"*-* obey equation (15) as before. T{u) obeys the following linear i?-matrix relation, 



{Taif3A'^),Ta2P2{'^)} = (^a2/3i W - ^"2/31 (^')) ^"1/32 + {T^iP2{v) - 71.1/32 (m)) 5a2P^) ■ 

Once again the spectral invariants contain the integrals of motion. Sklyanin's conjecture for the 
Gaudin model is. 

Conjecture 2 Let A and B he the same functions on GL{N) as above. Then the variables Xj andpj 
defined by the equations, 

B{T{x,)) = 0, = A{T{x,)) (69) 

have canonical Poisson brackets and, besides, are bound to the Hamiltonians by the relation 

det{pi-T{xi)) = {) . (70) 

This is true generically. To prove this one only needs to show that Xi and pj have the correct 
commutation relations. The proof works the same way as the quadratic case and will not be repeated 
(it can be found in Q). 

5.2 Other systems 

In this paper the procedure is illustrated only in the simplest cases of the magnetic chain and the 
Gaudin model, however the method is much more widely applicable. Since the equation M{u, C)Vo = 
and the calculation used to obtain A and B is independent of the particular i2-matrix algebra, the 



separation variables thus obtained satisfy condition (1) of section (1.3) (i.e. give separated equations) 
automatically. 

However it must be checked that they have the correct commutation relations and at present I do 
not have a general proof of this. Nevertheless it can be proved (by direct calculation) that one obtains 
the correct commutation relations for systems with i?-matrices of the following form 

R = ^fij{u)eij®eji (71) 

provided that the fij{u) = ^ + 0{l) as li — > 0. With the Cij as in equation ( pO| ) (which is clearly of this 
form). This allows the result to be extended from rational to trigonometric and elliptic /^-matrices. 
Applications in these systems and systems with dynamical /^-matrices will be discussed in a future 
publication. 
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6 Conclusion 



In this paper I have shown that genericahy the following polynomials in the T matrix coefficients 
maybe used to obtain separation variables as illustrated in the case of the magnetic chain and Gaudin 
model, 



n2 T->7V-2 ^»JV-il 

'det(M) 

a rpN-2 rp{N-l) 



A{T{u)) = (-iii2---iN-iTi^NTf^N---T!l_lN-J^(ij\ C^^) 



B{T{u)) = e,,,,...,^_^T,,^T^^---T:^^^T^;^:^ (73) 

in the sense that separation variables, with canonical Poisson brackets, are given by B{xi) = 0, 
gPi — p. — ^(xj) for the magnetic chain and by B{xi) = 0, pi = A{xi) for the Gaudin model. 

With A and B in this form it is clear to see that B has an SL(N — 1) symmetry corresponding 
to similarity transform which leave the last row and column fixed, and that A has an SL{N — 2) 
corresponding to those transformations which leave the first and last rows and columns fixed. Clearly 
in the construction and proof it does not matter which column we choose to call the first and last, 
and these associated ^'s and B^s also give rise to separation variables. The meaning of the symmetry 
and non-uniqueness is not yet clear. 

As remarked, B only gives the correct number of separation coordinates in the generic case. For 
non-generic cases the equation 

M{u,C)Vo = (74) 

must be used, it is desirable to extend the FBA to these cases. One would need to provide a method 
that told one how to decouple equation ( |4^ ) into an equation(s) for the Xi and an equation(s) for the 
PiS (in terms of the xis). The reducible systems of Kalnins et al.[ll, ^] provide examples of such 



problems and the author conjectures that the notion of irreducibility in these systems is related to 
this problem of non-degeneracy. 

The classical Functional Bethe Ansatz in its 2 x 2 version has been successfully applied to systems 
with a wide variety of i?-matrices including dynamical i?-matrices, as has its quantum counterpart. It 
is hoped that the SL{N) FBA proposed here and its quantum counterpart might prove equally useful. 
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